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PROJECT  OVERVIEW 


OBJECTIVES 

The  general  objective  of  the  project  is  fundamental  mathematical  modeling  of  a  complex  TSA 
system  with  electrothermal  desorption  step,  with  adsorbers  assembled  of  one  or  more 
cartridge-type,  radial  flow  fixed  beds  and  with  in-vessel  condensation  of  the  desorbed  vapor. 

During  the  first  stage  of  the  project,  a  1-D  mathematical  model  of  a  single  cartridge  without 
condensation  was  developed  in  Matlab  and  used  for  simulation.  During  the  second  stage  of 
the  project,  the  complete  TSA  system,  for  adsorbers  with  one  and  two  cartridges,  was 
modeled,  using  a  combination  of  Femlab  and  Matlab  software.  2-D  Femlab  models  that 
describe  adsorption,  electrothermal  desorption  and  electrothermal  desorption  accompanied 
with  condensation  of  the  desorbed  vapor  were  built  and  integrated  using  Matlab,  in  order  to 
describe  and  simulate  complete  TSA  cycles. 

The  objective  of  this,  third  phase  of  the  project,  is  to  develop  mathematical  models  of  another 
complex  configuration,  with  four  parallel  annular,  radial-flow  cartridge-type  adsorbers.  The 
analysis  of  this  configuration  shows  that  3-D  models  have  to  be  used  in  order  to  describe  it. 
The  models  for  adsorption,  electrothermal  desorption  and  electrothermal  desorption 
accompanied  with  condensation  of  the  desorbed  vapor  are  developed  again,  by  using 
COMSOL  Multiphysics  (Femlab)  3-D  modeling  tools. 


STATUS  OF  EFFORT 

Three  COMSOL  Multiphysics  3-D  models  were  built,  corresponding  the  three  stages  of  the 
TSA  process  with  electrothermal  desorption  taking  place  in  a  4-cartridge  adsorber: 
adsorption,  desorption  without  condensation  and  desorption  accompanied  with  condensation 
at  the  adsorber  wall.  The  models  incorporate  seven  application  modes  defined  in  a  number  of 
subdomains.  This  results  with  high  demands  regarding  the  computer  memory  and  computer 
time  needed  for  their  solution.  In  order  to  reduce  these  demands,  the  models  used  for 
simulations  were  defined  for  one  quarter  of  the  4-cartridge  adsorber,  assuming  that  all  four 
cartridges  are  identical.  Also,  the  solution  of  our  models  did  not  converge  directly.  A  three- 
step  iterative  solution  procedure  was  developed,  leading  to  the  solution  of  the  complete  3-D 
models.  Following  this  procedure  the  three  stages  of  the  TSA  process  were  simulated  for  a 
chosen  set  of  parameters. 


ACCOMPLISHMENTS  /NEW  FINDINGS 

The  accomplishments  of  this  project  can  be  perceived  from  two  aspects: 

1 .  The  developed  models  represent  a  mathematical  tool  which  can  be  used  for  simulation 
of  the  investigated  TSA  system,  its  analysis,  predicting  of  the  concentration  and 
temperature  profiles  in  the  system  and  its  optimization.  This  was  the  main  objective  of 
this  project. 

2.  The  developed  solution  procedure  in  COMSOL  Multiphysics  can  be  used  for  solving 
of  other  complex  3-D  models. 
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DETAILED  REPORT 


1.  Introduction 


The  idea  about  regeneration  of  used  adsorbents  by  direct  heating  of  the  adsorbent  particles  by 
passing  electric  current  through  them  (Joule  effect),  was  first  published  in  the  1970s  [1], 
Desorption  process  based  on  this  principle  was  later  named  electrothermal  desorption  [2],  It 
was  recognized  as  a  prospective  way  to  perform  desorption  steps  of  TSA  cycles.  At  the  same 
time,  fibrous  activated  carbon  was  recognized  as  a  very  convenient  adsorbent  form  for  its 
realization.  Electrothermal  desorption  has  some  advantages  over  conventional  methods, 
regarding  adsorption  kinetics  and  dynamics  [3,  4]  and  energy  efficiency  [2,  5].  Some 
industrial  applications  of  electrothermal  desorption  have  been  reported  recently  [6,  7], 
Nevertheless,  the  development  of  mathematical  description  of  processes  based  on 
electrothermal  desorption  hasn’t  been  following  the  development  of  the  process  itself,  so  far. 

A  new  TSA  process  with  electrothermal  desorption  step,  based  on  adsorbers  assembled  of  one 
or  more  annular,  cartridge-type,  fixed-beds,  with  in-vessel  condensation,  has  been  presented 
recently  [5,  8,  9].  Its  schematic  representation  is  given  in  Figure  1. 
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Figure  1.  Overall  schematic  representation  of  the  ACFC  adsorption  -  rapid  electrothermal  desorption  system 

(from  Ref.  [5]) 


The  aim  of  this  project  is  to  develop  rigorous  mathematical  models  of  this  system  for  different 
adsorber  configurations.  These  models  could  be  used  for  simulation,  analysis  and 
optimization  of  the  TSA  process  developed  by  Sullivan  [5], 


1.1.  Description  of  the  system 

The  adsorbers  of  the  TSA  system  shown  in  Figure  1  are  composed  of  one  or  more  annular, 
radial-flow  cartridge-type  adsorbent  beds.  A  single  cartridge  is  schematically  shown  in  Figure 
2.  It  is  formed  as  a  cylindrical  roll  of  activated  carbon  fiber  cloth  (ACFC),  spirally  coiled 
around  a  porous  central  pipe.  The  gas  flow  through  the  adsorber  is  in  the  radial  direction. 
During  the  desorption  step,  electric  current  is  passed  through  the  activated  carbon  cloth  in  the 
axial  direction,  causing  heat  generation,  heating  of  the  adsorbent  and  desorption. 
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The  cartridges  are  situated  in  a  cylindrical  adsorber 
vessel.  During  desorption,  the  hot  gas  steam  rich  with 
the  desorbed  vapor  comes  into  contact  with  the  cold 
vessel  wall  and  the  vapor  condenses.  In  that  way,  the 
adsorber  vessel  serves  as  a  passive  condenser.  The 
bottom  outlet  of  the  adsorber  vessel  is  funnel-shaped,  in 
order  to  facilitate  draining  of  the  condensed  liquid. 

The  structure  of  the  adsorbers  with  one,  two  and  four 
cartridges  is  shown  in  Figure  3.  In  principle,  the  flow 
directions  during  the  adsorption  and  desorption  steps 
can  be  different,  and  the  flow-rates  during  desorption 
are  considerably  lower  than  during  adsorption.  The 
flow  directions  shown  in  Figure  3  correspond  to  the 
desorption  step. 


Figure  2.  Schematic  representation  of 
the  annular  -  radial  flow  -  ACFC 

arknrhpnt  FpH 


Figure  3.  Adsorbers  with  one  (A),  two  (B)  and  four  (C)  annular,  radial-flow  cartridges 


1.2.  Brief  survey  of  the  results  of  the  first  two  phases  of  the  project 

Our  efforts  during  the  first  phase  of  the  project  were  focused  on  modeling  of  a  single¬ 
cartridge  adsorber,  without  condensation.  Two  deterministic  models  were  postulated:  the  first 
one  assuming  uniform  adsorbent  density  throughout  the  adsorbent  bed,  and  the  second  one 
taking  into  account  the  structure  of  the  bed  made  of  layers  of  ACFC.  Both  models  were 
obtained  as  complex  sets  of  coupled  nonlinear  PDEs,  ODEs  and  algebraic  equations. 
Numerical  procedures  for  solving  the  model  equations  were  established  and  used  for 
simulation  of  adsorption,  electroresistive  heating,  electrothermal  desorption  and  consecutive 
adsorption-desorption,  using  Matlab  software.  Similar  results  were  obtained  with  the  two 
models.  The  results  of  this  phase  are  presented  in  the  Final  Performance  Report  for  project 
FA8655-03-1-3010  [10]. 
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In  the  second  phase  of  the  project,  models  for  two  types  of  adsorbers,  one  with  only  one 
(Figure  3A),  and  the  other  with  two  cartridges  (Figure  3B),  have  been  developed.  For  each 
adsorber  type,  three  models  were  built,  in  order  to  describe  three  stages  of  a  compete  ESA 
cycle:  adsorption,  electrothermal  desorption  and  electrothermal  desorption  with  in- vessel 
condensation.  These  models  were  built  using  Femlab,  a  specialized  software  tool  for 
modeling  of  complex  systems  with  complex  geometry.  The  models  were  defined  in  2-D 
axially  symmetrical  geometry,  with  seven  application  modes.  In  order  to  describe  the 
complete  TSA  cycles,  the  models  for  the  three  stages  of  the  TSA  process  were  integrated,  by 
using  a  combination  of  Femlab  and  Matlab.  The  models  were  successfully  used  for  simulation 
of  separate  stages  of  the  process  and  of  the  complete  TSA  cycles,  as  well  as  for  their 
optimization.  The  results  of  this  phase  are  presented  in  the  Final  Performance  Report  for 
project  FA8655-04-1-3053  [11], 


1.3.  Outline  of  the  third  phase  of  the  project 

The  third  phase  of  the  project  was  focused  on  mathematical  modeling  of  the  four-cartridge 
adsorber,  presented  in  Figure  3C.  This  configuration  has  been  described  in  detail  in  Refs.  10 
and  12.  This  modeling  is  performed  using  COMSOL  Multiphysics  software  (previously 
Femlab).  Owing  to  the  complex  configuration  of  this  adsorber  design,  3-D  models  have  to  be 
used,  leading  to  a  number  of  issues  regarding  computer  memory  requirements  and 
convergence  problems. 

Similarly  to  the  2-D  models  developed  in  the  second  phase,  three  3-D  models  are  needed,  to 
describe  three  steps  of  the  complete  cycle: 

adsorption, 

desorption  without  condensation, 
desorption  accompanied  with  condensation. 

The  following  application  modes  were  incorporated  in  these  models:  Non-isothermal  flow, 
Brinkman  Equations,  Convection  and  Diffusion,  Diffusion,  Convection  and  Conduction,  Heat 
transfer  by  Conduction  and  Conductive  Media  DC. 

Solving  of  these  very  complex  3-D  models  is  a  rather  demanding  and  difficult  task.  It  assumes 
optimization  of  the  model  structure  and  geometry  definition,  the  generated  mesh  and  the 
choice  of  the  solver,  in  regard  to  the  used  hardware,  computer  time  and  convergence  of  the 
solution. 

All  these  issues  were  attended  in  this  phase  of  the  project. 


2.  Modeling  of  the  four-cartridge  adsorber,  using  COMSOL  Multiphysics 

Three  3-D  models  were  built,  for  the  three  steps  of  the  complete  TSA  cycle: 

•  Model_A3  -  for  adsorption 

•  Model_D3  -  for  desorption  without  condensation 

•  Model_DC3  -  for  desorption  with  condensation 
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Model  assumptions: 

The  same  assumptions  as  for  the  2-D  models  derived  in  the  second  phase  of  the  project  [11] 
were  used  in  setting  up  the  3-D  models: 

•  The  adsorbent  beds  are  treated  as  homogeneous,  with  uniform  adsorbent  porosity  and 
density. 

•  The  mass  and  heat  transfer  resistances  on  the  particle  scale  are  neglected. 

•  The  fluid  phase  is  treated  as  an  ideal  gas  mixture  of  the  inert  and  the  adsorbate. 

•  All  physical  parameters  and  coefficients  are  considered  as  constants. 

•  The  electric  resistivity  of  the  ACFC  adsorbent  is  temperature  dependent.  Linear 
temperature  dependence  is  assumed,  based  on  experimental  results  reported  in  Ref. 
[5], 

•  The  electric  power  during  electrothermal  desorption  is  supplied  under  constant  voltage 
conditions. 

•  The  condensation  at  the  adsorber  wall  is  dropwise.  This  assumption  is  based  on 
experimental  observations  reported  in  Ref.  [5], 

•  The  volume  of  the  condensed  liquid  is  neglected,  i.e.,  it  is  assumed  that  the  liquid 
drops  don’t  influence  the  gas  flow. 

•  The  heat  resistance  and  heat  capacity  of  the  adsorber  wall  are  neglected,  so  that  the 
wall  temperature  is  equal  to  the  temperature  of  the  environment. 

•  Initially,  the  adsorbate  concentrations  in  both  phases  are  in  equilibrium  and  uniform 
throughout  the  adsorbent  bed.  The  temperatures  of  both  phases  are  initially  equal  and 
uniform  throughout  the  adsorbent  bed. 


2.1.  Definition  of  the  model  geometry 

In  the  second  phase  of  the  project  we  built  Femlab  models  for  adsorbers  with  one  and  two 
cartridges  and  we  used  2-D  axially  symmetrical  space  [11],  Nevertheless,  the  four-cartridge 
adsorber  shown  in  Figure  3C  is  not  axially  symmetrical,  so  it  has  to  be  modeled  in  3-D  space. 
In  COMSOL  Multiphysics,  the  only  option  for  3-D  modeling  is  to  use  the  Descartes 
coordinate  system  (x,  y,  z). 

The  gas-flow  and  electric  current  directions  in  the  four-cartridge  adsorber,  as  defined  in  Ref. 
[12],  are  shown  in  Figure  4.  This  figure  shows  that  the  gas-flow  directions  differ  during 
adsorption  and  desorption.  During  adsorption,  two-by-two  cartridges  work  in  parallel:  one 
pair  of  cartridges  serves  for  introducing  the  feed  stream,  and  the  other  for  removal  of  the  gas 
stream  from  the  adsorber.  The  two  pairs  of  cartridges  are  connected  in  series.  From  the 
modeling  point  of  view,  the  cartridges  used  for  introducing  the  gas  stream  into  the  adsorber 
are  identical,  as  well  as  the  cartridges  used  for  withdrawing  the  gas  stream  from  the  adsorber, 
but  the  two  pairs  of  cartridges  are  different.  On  the  other  hand,  during  desorption,  all  four 
cartridges  work  in  parallel,  with  gas  inlet  from  the  top  ends  of  the  central  tubes  and  gas  outlet 
through  the  bottom  end  of  the  outer  tube  in  which  all  four  cartridges  are  situated  (the  adsorber 
vessel).  The  condensed  liquid  during  desorption  leaves  the  adsorber  through  the  same  outlet. 
From  the  modeling  point  of  view,  all  four  cartridges  are  identical,  regarding  momentum,  mass 
and  heat  transfer.  The  electric  current  flows  are  also  shown  in  Figure  4.  For  the  design  shown 
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in  this  figure,  all  four  cartridges  are  connected  in  series.  From  the  modeling  point  of  view,  all 
four  cartridges  are  different,  regarding  the  electric  current  flow. 


Vessel  lid  is  constructed  of 
non-conducting  material. 


Gas  out 


Nitrogen 


Gas  in 


Gas  flows  through  the 
cartridges  two  in  series 
and  two  in  parallel. 


Temperature 
measurement  locations 

1.  ACFC 

2.  Exit 


The  bottoms  of  the  cartridges  are 
sealed  with  a  conducting  plug. 
Adjacent  cartridges  are  electrically 
bridged.  The  left- most  cartridges  are 
electrically  bridged  at  the  top  (not 
shown). 


Condensate  and  N2 


Electric  current  flows  through  all 
cartridges  in  series. 


ADSORPTION 


REGENERATION 


Figure  4.  Four-cartridge  adsorber  [12] 


In  principle,  the  complete  adsorber  geometry  with  four  cartridges  should  be  used,  in  order  to 
model  the  adsorber  in  which  all  four  cartridges  are  different  (e.g.  because  of  different  electric 
current  flows  through  them).  Nevertheless,  3-D  models  are  very  complex  and  very  demanding 
regarding  computer  memory  and  computer  time.  For  that  reason,  we  explored  two  other 
options:  to  model  a  four-cartridge  adsorber  in  which  two  cartridges  can  be  considered  as 
equal,  and  a  four-cartridge  adsorber  in  which  all  four  cartridges  are  equal. 


2.1.1.  A  3-D  model  of  the  whole  adsorber 

The  most  general  and  most  complex  geometry  is  the  one  representing  the  whole  adsorber  with 
4  cartridges.  It  is  shown  in  Figure  5,  in  3-D  perspective  and  in  the  2-D  x-y  plane.  The 
complete  adsorber  consists  of  9  sub-domains  (shown  in  the  2-D  view  in  Fig.  5): 

Sub-domains  EXT1,  EXT2,  EXT3  and  EXT4  correspond  to  the  adsorbent  beds  of 
the  4  cartridges; 

Sub-domain  EXT5  corresponds  to  the  adsorber  space  around  the  4  cartridges; 

Sub-domains  EXT6,  EXT7,  EXT8  and  EXT9  correspond  to  the  central  tubes  of 
the  4  cartridges. 
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Figure  5.  Geometry  definition  for  the  whole  four-cartridge  adsorber:  3-D  view  (left)  and  2-D  x-y  view  (right) 

If  “coarser  mesh”  is  chosen  as  a  mesh  parameter,  with  all  other  default  parameter  settings,  the 
mesh  corresponding  to  this  geometry  contains  24481  3-D  elements. 


2.1.2.  A  3-D  model  of  one  half  of  the  adsorber 

If  the  adsorber  contains  two  pairs  of  identical  cartridges,  the  model  geometry  could  be 
simplified  by  defining  only  one  half  of  the  adsorber  and  assuming  that  the  other  half  is 
completely  symmetrical  to  it.  This  geometry  is  shown  in  Figure  6. 


Figure  6.  Geometry  definition  for  one  half  of  the  four-cartridge  adsorber:  3-D  view  (left)  and  2-D  x-y  view 
(right) 


For  this  geometry,  the  number  of  sub-domains  reduces  to  5  (shown  in  the  2-D  view  in  Fig.  6): 

Sub-domains  EXT1  and  EXT2  correspond  to  the  adsorbent  beds  of  the  two 
cartridges; 

Sub-domain  EXT3  corresponds  to  the  space  around  the  cartridges; 

Sub-domains  EXT4  and  EXT5  correspond  to  the  central  tubes  of  the  two 
cartridges. 

For  “coarser  mesh”  setting,  the  mesh  for  this  geometry  consists  of  12118  3-D  elements. 


2.1.3.  A  3-D  model  of  one  quarter  of  the  adsorber 

If  all  four  cartridges  can  be  considered  as  identical,  the  model  geometry  can  be  further 
simplified,  by  defining  only  one  quarter  of  the  adsorber  and  assuming  symmetry  both  in  the  x- 
and  in  the  y-direction.  This  geometry  is  shown  in  Figure  7  (the  main  dimensions  of  the 
adsorber  are  also  shown). 


Figure  7.  Geometry  definition  for  one  quarter  of  the  four-cartridge  adsorber:  3-D  view  (left)  and  2-D  x-y 
view  (right) 


In  this  case,  three  sub-domains  have  to  be  defined  (show  in  the  2-D  view  in  Fig.  7): 

Sub-domain  EXT1  -  the  adsorbent  bed; 

Sub-domain  EXT2  -  the  space  around  the  cartridges; 

Sub-domain  EXT3  -  the  central  tube  of  the  cartridge. 

The  mesh  for  this  geometry,  for  “coarser  mesh”  setting,  contains  6126  elements. 

All  three  geometries  shown,  in  Figures  5,  6  and  7,  were  defined  for  a  laboratory  scale 
adsorber  of  the  following  dimensions: 
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-  H=  10  cm  (height  of  the  adsorber) 
ri=0.95  cm  (radius  of  the  central  tubes) 

ri=  1 .5  cm  (outer  radius  of  the  adsorbent  beds) 
r3= 5  cm  (radius  of  the  adsorber  vessel) 

-  xo=Vo=2  cm  (position  of  the  cartridge  center  in  the  x  and  y  directions) 


2.2.  The  application  modes  used 

For  modeling  of  the  four-cartridge  adsorber,  the  same  built-in  application  modes  are  used,  as 
for  the  one-  or  two-cartridge  adsorbers  (Ref.  11).  For  the  most  complex  process  - 
electrothermal  desorption,  in  which  simultaneous  gas  flow  throw  tubes  and  adsorbent  beds, 
mass  and  heat  transfer  and  heat  generation  by  the  Joule  effect  are  coupled,  seven  application 
modes  have  to  be  incorporated  in  the  COMSOL  Multiphysics  model: 

1 .  N on-isothermal  flow  -  for  defining  the  momentum  balances  in  the  central  tubes  and 
outer  space; 

2.  Brinkman  equations  -  for  defining  the  momentum  balances  for  the  adsorbent  beds; 

3.  Convection  and  Conduction  -  for  defining  the  heat  balances  for  the  gas  phase; 

4.  Heat  transfer  by  Conduction  -  for  defining  the  heat  balances  for  the  solid  phase; 

5.  Conductive  Media  DC  -  for  defining  heat  generation  by  Joule  effect; 

6.  Convection  and  Diffusion  -  for  defining  the  mass  balances  for  the  gas  phase; 

7.  Diffusion  -  for  defining  the  mass  balances  for  the  solid  phase. 

The  application  modes  are  active  in  the  sub-domains  according  to  this  scheme: 

•  Convection  and  Conduction  and  Convection  and  Diffusion  are  active  in  all  sub- 

domains; 

•  N on-isothermal  flow  is  active  in  sub-domains  corresponding  to  the  central  tubes  of 
the  cartridges  and  the  space  around  the  cartridges; 

•  Brinkman  equations,  Heat  transfer  by  Conduction,  Conductive  Media  DC  and 
Diffusion  are  active  only  in  the  sub-domains  corresponding  to  the  adsorbent  beds. 


2.3  The  underlying  equations 

In  principle,  the  phenomena  taking  place  in  the  four-cartridge  TSA  system  are  the  same  as  in 
the  one-  and  two-cartridge  systems.  Nevertheless,  as  a  result  of  using  3-D  Descartes  instead  of 
2-D  radial  coordinate  system,  the  model  equations  are  not  exactly  the  same.  Although  the 
model  equations  and  their  boundary  conditions  for  different  steps  of  the  TSA  process  are 
somewhat  different,  as  electrothermal  desorption  involves  all  phenomena  and  the 
corresponding  application  modes,  here  is  an  overview  of  the  PDEs  underlying  the  FEMLAB 
model  for  electrothermal  desorption: 
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1)  Momentum  balances  and  continuity  equations  for  the  inert  gas: 

For  the  central  tubes  and  the  space  around  the  cartridges'. 
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For  the  adsorbent  beds : 
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The  following  notations  are  used  in  these  equations:  t  -  time,  x,  y,  z  -  space  coordinates,  u,  v 
and  w  -  gas  velocities  in  the  x,  y  and  z  direction,  respectively,  p  -  pressure,  pg  -  inert  gas 
density,  k  -  adsorbent  bed  permeability  and  p  -  dynamic  viscosity  of  the  carrier  gas  (inert). 


2)  Adsorbate  balance  for  the  gas  phase: 

For  the  central  tubes : 
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For  the  adsorbent  beds : 


3(P»C)  .  J  „  ,  <5(p,C)  „  ,  d(psC)  „  ,  a(p,C) 


dt 


■+v 


■D 


dx 


D 


D 


mz  \b 


dz 


=  kma(C*-C) 


d(psC )  d(pgC)  d(/?gC) 
s  +w- 


dx 


-\-  U  - 


dy 


dz 


(10) 


For  the  space  around  the  cartridges'. 


d(PgC)  J 
-  + V 


dt 


V 


D  I  W 

mxko  dx 


D  I  ^£)  D  I 

dy  wz|i-°  dz 


e(Pgc)  d(Pgc)  d(Pgc) 

u - - —  +  u - - —  +  w  s 


V 


dx 


dy 


dz 


J 


(11) 


In  these  equations  C  is  adsorbate  concentration  in  the  gas  phase,  defined  as  molar  ratio,  i.e.,  in 
moles  of  adsorbate  per  mole  of  inert  gas  and  C*  the  gas  concentration  in  equilibrium  with  the 
solid  phase.  Dmx,  Dmy  and  Dmz  are  mass  dispersion  coefficients  in  the  x,  y  and  z  direction, 
respectively,  km  mass  transfer  coefficient  and  a  specific  surface  area  of  the  adsorbent  for  mass 
and  heat  transfer.  The  subscript  i,o  corresponds  to  the  inner  (central)  tubes  and  outer  space 
and  b  to  the  adsorbent  beds. 


3)  Heat  balances  for  the  gas  phase: 

For  the  central  tubes  and  the  space  around  the  cartridges : 
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For  the  adsorbent  beds : 
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+  \a(T,-T) 


In  equations  (12)  and  (13)  Tg  and  Ts  are  the  gas  and  solid  temperature,  respectively, 
D'u  ,D^  and  D^8  heat  diffusivities  of  the  gas  phase,  in  the  x,  y  and  z  direction,  respectively,  hb 

heat  transfer  coefficient  within  the  adsorbent  bed  and  cpg  and  cpv  molar  heat  capacities  of  the 
inert  gas  and  adsorbate  vapor,  respectively. 


4)  Electric  current  balance  for  resistive  heating: 

J  dU  1  dU  1  dU\  A 
-V - + - + - =0 

vp  dx  p  dy  p  dz  J 


(14) 
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In  this  equation  U  is  electric  potential  and  p  electric  resistivity  of  the  adsorbent,  which  is 
temperature  dependent.  Linear  temperature  dependence  of  the  electric  resistivity,  which  was 
obtained  experimentally  in  Ref.  5,  was  used  in  our  models: 

p  =  Po(i  +  6(7;-rfi))  (15) 


5)  Adsorbate  balance  for  the  solid  phase  within  the  adsorbent  bed: 

p  M  =  kma(C-C)  (16) 

dt 

where  q  is  the  adsorbate  concentration  in  the  solid  phase  and  p b  the  solid  phase  density.  The 
Dubinin-Radushkevich  equation,  valid  for  adsorption  of  hydrocarbons  on  activated  carbon  [5] 
is  used  to  describe  the  adsorption  equilibrium.  This  equation  can  be  put  into  the  following 
form: 


C  = 


p  expl 


-In 


'MAq' 

Pa% 


(17) 


where  Wo  is  the  volume  of  the  micropores  per  unit  mass  of  adsorbent,  E-  the  adsorbate  energy 
of  adsorption  and  Ma  and  pj  the  adsorbate  molar  mass  and  density,  respectively.  p°  is  the  the 
adsorbate  saturation  pressure.  In  our  models,  it  was  calculated  using  the  Wagner  equation: 


r  VPAX  +  VPBx1'5  +  VPcx3  +  VPDx  ' 

[PcJ 

l  J 

where  x  =  \-(TjTc) ,  pc  and  Tc  are  the  critical  pressure  and  temperature  and  VP  a,  VPb,  VPc 
and  VP/)  are  the  Wagner  constants. 


6)  Heat  balance  for  the  solid  phase  within  the  adsorbent  bed: 
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where  (&Qe,/dV)  is  the  electric  power  supply  per  unite  volume  of  the  adsorbent  bed: 


8 Ol  =  I 

dV  p 


v  dx  j 


+ 
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In  the  case  of  adsorption  (SQel/dV)  =  0.  In  equation  (19)  (AHads)  is  the  molar  heat  of 

adsorption,  D1"  ,D^  and  the  heat  diffusivities  of  the  solid  phase  in  the  x,  y  and  z  direction, 

respectively,  and  cps  and  cpi  the  heat  capacities  of  the  solid  and  of  the  liquid  adsorbate, 
respectively. 
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Boundary  conditions: 

The  boundary  conditions  listed  bellow  correspond  to  the  gas  flow  dirrections  defined  for  the 
desorption  step  in  Figure  4,  and  for  one  quarter  of  the  four-cartridge  adsorber. 

-  For  the  bottom  surface  of  the  central  tube: 

z  =  0,  (x-x0)2  +  (y  —  y0)2  <  q2 :  u=  0,  v  =  0,  w  =  0, 


■D? 


—  =  0,  -  —  =  0 

,o  Qz  mz"’°  Qz 

-  For  the  top  surface  of  the  central  tube: 
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-  For  the  cylindrical  surface  of  the  central  tube: 
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-  For  the  bottom  surface  of  the  adsorbent  bed: 
z  =  0,  rj2  <  (x-x0)2  +(y  —  y0)2  <  r22  :  u  =  0,  v  =  0,  w  =  0, 
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(24) 


-  For  the  top  surface  of  the  adsorbent  bed: 

Z  =  H,  q2  <(x-x0)2  +(y-y0Y  <r2  :  u  =  0,  v  =  0,  w  =  0, 
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-  For  the  outer  cylindrical  surface  of  the  adsorbent  bed: 

(x  -  x0)2  +  (y  -  y  0)2  -r2 ,  ze(0,H):  p\b=p\ot,  u\b=u\ot,  v\b=v\ot,  4  =  H 
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-  For  the  bottom  surface  of  the  space  around  the  cartridges: 
z  =  0,  (x  -  x0)2  +  (x  -  y0 )2  >  r2  n  x2  +  y2  <  r2  n  x  e  (0 ,r3),  y  e  (0,r3) : 
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-  For  the  top  surface  of  the  space  around  the  cartridges: 
z  =  H,  (x  -  x0)2  +  (x  -  y0)2  >  r2  n  x2  +  y2  <  r2  n  x  e  (0  ,r3),  y  e  (0,r3)  : 

u=  0,  v  =  0,  w  =  0,  ~P>hf 
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-  For  the  adsorber  wall: 

x2+y2=r2,  ze(0,H ):  w  =  0,  v  =  0,  w  =  0, 
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(29) 


For  the  case  of  desorption  accompanied  with  condensation,  the  boundary  condition  (29)  has  to 
be  changed  in  the  following  way  [11]: 


x2+y2=r2,  ze(0,H ):  u  =  0,  v  =  0,  w  =  0, 
C  =  Csat(TJ,  ~d: 
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In  equations  (21-29),  r\,  rp  r 3,  xo,  yo  and  H  are  the  adsorber  dimensions  which  have  been 
defined  in  Figure  7.  G  is  the  molar  flow-rate  of  the  gas  inert,  Co  -  the  supply  voltage  (0  for 
adsorption)  and  Sb  -  bed  porosity.  The  subscript  in  denotes  inlet,  it  the  inner  (central)  tube,  ot 
the  space  around  the  cartridges,  a  the  ambient  and  w  the  adsorber  wall.  (A Hcon(j)  is  the  molar 
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heat  of  condensation  and  Csat(Tw)  the  saturation  concentration  corresponding  to  the  wall 
temperature  Tw  ,  which  was  calculated  using  the  Wagner  equation. 

Although  the  heat  and  mass  transfer  coefficients,  dispersion  coefficients,  etc.,  can  be  defined 
as  functions  of  the  dependent  variables  (velocities,  concentrations,  temperatures),  in  our 
current  models  constant  values  were  used  in  order  to  reduce  the  convergence  problems.  The 
adsorption  isotherm  relation  (q=f(C,  Ts))  and  the  temperature  dependence  of  the  electric 
resistivity  were  included  in  the  models  (equations  (17)  and  (15),  respectively). 


3.  Solution  of  the  models 

The  Models  of  the  four-cartridge  adsorber  were  developed  and  solved  using  a  personal 
computer  with  AMD  ATHLON™  64  3000+  processor  and  with  4  GB  RAM. 

Our  3-D  models  are  very  complex,  regarding  both  the  geometry  and  the  number  of 
application  modes  used.  Even  with  a  “coarser”  mesh,  the  number  of  degrees  of  freedom 
obtained  is  very  high.  This  number  was  so  high,  that  the  models  for  the  whole  and  for  one  half 
of  the  four-cartridge  adsorber  (Figs.  5  and  6)  could  not  be  solved  at  all,  owing  to  the  lack  of 
computer  memory.  Because  of  that,  we  focused  our  efforts  on  solving  the  model  of  a  quarter 
of  the  adsorber,  assuming  model  symmetry  in  the  x  and  y  directions  (Fig.  7). 


3.1.  The  mesh  parameters 

In  order  to  reduce  the  memory  requirements,  the  following  mesh  parameters  were  used: 
Predefined  mesh  sizes:  “coarser” 

Maximum  element  size  scaling  factor:  1 .9 
Element  growth  rate:  1 .7 

With  these  mesh  parameters  the  following  mesh  statistics  was  obtained: 

Number  of  mesh  elements:  2237 
Minimum  element  quality:  0.2578. 

The  mesh  with  these  characteristics  is  presented  in  Figure  8. 
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Figure  8.  The  mesh  for  one  quarter  of  the  four-cartridge  adsorber:  3-D  view  (left)  and  2-D  x-y  view  (right) 


3.2.  The  solver 

The  three  models  were  solved  using  the  direct  SPOOLS  solver,  with  default  solver 
parameters.  The  following  values  of  the  absolute  and  relative  tolerance  and  the  time  step  were 
used: 

For  Model_A3  (adsorption): 

•  Absolute  tolerance:  0.001 

•  Relative  tolerance:  0.01 

•  Time  step:  20  s 

For  Model_D3  and  Model_DC3  (desorption  without  and  with  condensation): 

•  Absolute  tolerance:  0.01 

•  Relative  tolerance:  0.1 

•  Time  step:  10  s 

The  solution  of  Model_D3  and  Model_DC3  did  not  converge  with  lower  tolerances  and 
longer  step  time. 


3.3.  The  solution  procedure 

Owing  to  their  great  complexity,  the  complete  models  couldn’t  be  solved  directly,  i.e.  the 
solution  usually  diverged.  The  main  difficulty  was  to  solve  the  nonisothermal  Navier-Stokes 
equations  simultaneously  with  the  heat  and  mass  balances,  and  especially  to  start  the  solution 
process  with  zero  initial  values  of  the  velocities.  In  order  to  overcome  these  problems,  a  step- 
by-step  solution  procedure  was  developed: 
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Step  1.  A  simplified  model,  in  which  only  Non-Isothermal  Flow  and  Brinkman  Equations 
application  modes  are  active,  is  solved  using  the  stationary  nonlinear  solver.  In  this  way  the 
steady  state  values  of  the  velocities  and  pressures  in  the  system  are  calculated,  for  the  case  of 
inert  isothermal  flow.  In  reality,  the  steady-state  velocities  and  pressures  are  established  very 
fast,  during  the  first  3  seconds  [1 1],  so  using  a  stationary  solver  is  appropriate.  In  this  step,  the 
initial  values  of  the  solved  variables  (velocities  in  all  three  sub-domains  and  x,  y  and  z 
direction  and  pressures  in  all  three  sub-domains)  are  taken  from  the  Initial  value  expression. 
The  values  of  the  variables  not  solved  in  this  step  (temperatures  and  concentrations  of  both 
phases)  are  taken  from  Initial  value  frame.  The  solution  of  this  stage  is  stored.  The  number  of 
degrees  of  freedom  for  this  calculation  was  14599. 

Step  2.  The  model  is  solved  with  the  following  application  modes:  Convection  and 
Conduction,  Heat  Transfer  by  Conduction,  Convection  and  Diffusion,  Diffusion  and 
Conductive  media  DC  (this  mode  is  not  used  in  Model_A3).  The  application  modes  Non- 
Isothermal  Flow  and  Brinkman  Equations  are  inactive.  In  this  step  the  heat  and  mass  balances 
are  solved,  using  the  time  depended  solver,  with  fixed  velocities  and  pressures,  calculated 
previously  in  Step  1 .  The  initial  values  of  the  variables  solved  in  this  step  (temperatures  and 
concentrations  of  both  phases)  were  taken  from  the  Initial  value  frame  and  the  values  of  the 
variables  not  solved  in  this  step  (velocities  and  pressures)  from  the  Stored  solution.  At  the 
end,  the  solution  at  time  0  was  stored.  The  number  of  degrees  of  freedom  was  10508  for 
Model_A3  and  12425  for  Model_D3  and  Model_DC3. 

Step  3.  The  integral  model  is  solved  using  the  time  depended  solver,  with  all  application 
modes  active.  The  initial  values  of  all  variables  are  taken  from  the  Stored  solution  at  time  0. 
For  faster  and  more  stable  convergence  this  step  is  performed  twice.  The  first  run  is 
interrupted  and  the  solution  at  time  0  is  stored,  and  a  new  run  is  performed  with  the  initial 
values  of  all  variables  taken  from  that  stored  solution.  If  necessary,  larger  tolerance  values  can 
be  used  in  the  first  run.  The  number  of  degrees  of  freedom  in  this  step  was  25107  for 
Model_A3  and  27024  for  Model_D3. 

Step  3  for  Model_DC3  did  not  converge,  even  when  the  tolerance  was  further  increased. 


4.  Simulation  results 

The  simulations  were  performed  for  a  laboratory  scale  adsorber,  with  the  dimensions  defined 
in  Section  3.1. 

The  model  parameters  correspond  to  the  following  system: 

-  Adsorbate:  methyl  ethyl  ketone  (MEK) 

-  Adsorbent:  American  Kynol  ACC-5092-20  (activated  carbon  fiber  cloth  -  ACFC 
material) 

-  Carrier  gas:  nitrogen. 

The  bed  permeability  was  calculated  from  the  experimentally  measured  pressure  drop  of  the 
ACFC  bed  given  in  Ref.  5.  The  numerical  values  of  most  of  the  other  model  parameters  are 
also  based  on  Ref.  5.  A  complete  list  of  the  model  parameters  used  in  these  simulations  is 
given  in  Table  1. 
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Table  1.  The  model  parameters  used  for  simulation 


Parameter 

Notation  and  value  | 

Adsorption 

Desorption  ] 

Adsorber  dimensions 

r,=0.95  cm  r2=l  .5  cm 

r3= 5  cm  x0=2cm,  y0= 2cm 

H=10 cm 

Bed  porosity 

86=0.72 

Specific  surface  area 

a=  13.65  cm2/cm3 

Ambient  pressure 

Atm^lO  1-325  kPa 

Inlet  adsorbate  concentration  in  the  gas  phase 

Cin= 0.001  mol/mol 

C/w=0  mol/mol 

Initial  adsorbate  concentration  in  the  gas  phase 

Cp= 0  mol/mol 

C^=0.001  mol/mol 

Initial  adsorbate  concentration  in  the  solid  phase 

qp= 0  mol/g 

^=0.00465  mol/g 

Saturation  concentration 

C^=0.0147  mol/mol  | 

Electric  voltage  supply  U0 

U0= 0  V 

U0= 3V 

Inert  gas  flow-rate 

G=0.03  mol/s 

G=0.01  mol/s 

Inlet  gas  temperature 

7^=293.15  K 

Ambient  temperature 

7>293.15  K 

Wall  temperature 

7^=293.15  K 

Mass  transfer  coefficient  in  the  bed 

km=0.0003  mol/(cm2s) 

Bed  to  central  tube  mass  transfer  coefficient 

kml=0.0003  mol/(cm2s) 

Bed  to  annular  tube  mass  transfer  coefficient 

km2=0.0003  mol/(cm2s) 

Heat  transfer  coefficient  in  the  bed 

hb=0.0\  W/(cm2s) 

Bed  to  central  tube  heat  transfer  coefficient 

hsl= 0.02  W/(cm2s) 

Bed  to  central  tube  heat  transfer  coefficient 

hs2= 0.012  W/(cm2s) 

Heat  transfer  coefficient  defining  heat  losses 

hm= 5x1 0'5  W/(cm2K) 

Dispersion  coefficient  of  the  gas  phase  in  the  central  and  annular  tubes 

^Tnx,y„zl  it- ^mx,y,zl  ot— 1  mol/(cms) 

Heat  diffusivity  of  the  gas  phase  in  the  central  and  annular  tubes 

Ax.v.zhgl  it=Ax.v.zhgl  „,=1  W/(Kcm) 

Heat  diffusivity  of  the  solid  phase 

Afe=1.2xl0-3  W/(cmK) 

Molar  heat  of  adsorption 

(~AHads)=66200  J/mol 

Molar  heat  of  condensation 

(-A//cW)=31230  J/mol 

Adsorbent  bed  density 

p6=0.221  g/cm3 

Adsorbate  density 

p^=0.81  g/cm3  — — —  — 

Adsorbate  molar  mass 

Ma=72A01  g/mol 

Inert  gas  molar  mass  (nitrogen) 

Mb= 28.02  g/mol 

Specific  heat  capacity  of  liquid  adsorbate 

cd/=157.9  J/(molK) 

Specific  heat  capacity  of  the  inert  gas 

cv=29A3  J/(molK) 

Heat  capacity  of  the  solid  phase 

CpS=0.71  J/(gK) 

Specific  heat  capacity  of  gaseous  adsorbate 

CpV=  100.9  J/(molK) 

Dynamic  viscosity  of  the  inert  gas 

p=  1.78x1  O'4  g/cm/s 

Inert  gas  thermal  conductivity 

*=2.6x1  O'4  W/cm/K 

Electric  resistivity  at  referent  temperature  TR 

p0=0.202  Hem 

Temperature  coefficient  of  the  bed  electrical  resistivity 

*=-2xl  0'3  1/K 

Critical  temperature  of  the  inert  gas  (nitrogen) 

TcB=  126.2  K 

Critical  pressure  of  the  adsorbate 

pcA= 4.26x1 0+6  Pa 

Critical  pressure  of  the  inert  gas  (nitrogen) 

pcB= 3.39xl0+6  Pa 

Wagner  constants  for  the  adsorbate 

F^=-7.71476,  VPb=\.1\06\ 
FPc^-3.6877,  FPd=-0.75169 

Total  volume  of  micropores  (D-R  equation) 

JF0= 0.748  cm3/g 

Boltzmann’s  constant 

*b=1.38048xl0"23  J/K 

Adsorption  energy  of  the  adsorbate  (D-R  eq.) 

£'=14.43xl0+3  J/mol 
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COMSOL  Multiphysics  offers  a  large  number  of  options  for  graphical  representation  of  the 
simulated  results.  The  following  results  were  chosen  for  presentation  in  this  report: 

Velocity  field  in  the  adsorber; 

Gas  and  solid  concentration  fields  in 
the  adsorber; 

Gas  and  solid  temperature  fields  in 
the  adsorber. 

For  all  three  models,  these  results  are  shown  in 
the  form  of  3-D  slice  plots,  corresponding  to  3 
different  times.  The  arrows,  showing  the 
direction  of  the  gas  flow  are  also  shown  in  these 
3-D  plots.  Also,  2-D  plots  are  given, 
corresponding  to  two  cross-sections  of  the 
adsorber,  one  corresponding  to  x=y,  and  the 
other  to  z=H!2  (planes  A  and  B  in  Figure  9). 

Time  profiles,  in  1-D  diagrams,  corresponding 
to  5  points  in  the  adsorber,  also  defined  in  Figure  9,  are  also  given.  These  points  lie  at  the 
intersection  of  planes  A  and  B,  two  of  them  within  the  adsorbent  bed,  two  in  the  space  around 
the  cartridges  and  one  on  the  axes  of  the  central  tube. 

Two  additional  parameters,  very  important  for  the  electrothermal  desorption  process,  were 
calculated  and  presented:  the  amount  of  the  condensed  liquid  and  the  used  electric  energy. 


Figure  9.  Definition  of  the  planes  and  points 
used  for  2-D  and  1-D  presentation 


Calculation  of  the  amount  of  the  condensed  liquid: 

The  flux  of  the  condensed  liquid  is  calculated  from  the  gas  concentration  gradient  at  the 
adsorber  wall: 


Jcond{^y^Zd)  =  -Dm 


I 

1,0  dx 


D 


my 


8C_ 
-°  dy 


* +y  =? 


(30) 


It  should  be  noticed  that,  owing  to  the  3-D  geometry,  this  flux  is  a  function  of  all  three 
coordinates  and  time.  The  condensation  rate  is  calculated  by  integration  of  the  flux  over  the 
whole  surface  of  the  adsorber  wall,  and  is  a  function  of  time  only: 

kond  (0  =  4  Jj  Jcond  ( x,y,z ,  t)dxdydz  (3 1 ) 

x2+y2=d 
0  <z<H 


The  total  amount  of  the  condensed  liquid  is  obtained  by  integrating  the  condensation  rate  over 
time: 


4  ond=  \ko,Jt)dt  (32) 

0 

^cond  is  the  total  time  of  condensation. 
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Calculation  of  the  used  electric  energy: 

It  is  assumed  that  the  used  electric  energy  is  equal  to  the  Joule  heat  produced  in  the  adsorbent 
material.  This  energy  per  unit  volume  and  unit  time  is  defined  by  equation  (20). 

The  used  electric  power  per  cartridge  is  obtained  by  integration  over  the  adsorbent  bed 
volume: 


JJJ 


i 


n  <(*-*o)  +(y-y<>)  <ri 

0  <z<H 


A(1  +  b(Ts-TR)) 


f  dU} 

\dx  j 


2  f  arrV 


+ 


dU 


+ 


rdU} 

V  dz  j 


dxdydz 


(33) 


and  the  total  used  energy  during  desorption,  by  integration  over  time  and  multiplying  by  4: 


^ des 

Qel  =  |  Qel  C)dt 


(34) 


x des  is  the  total  desorption  time. 


4.1.  Simulation  of  adsorption 

The  simulation  of  adsorption  was  performed  using  Model_A3.  In  order  to  reduce  the  memory 
and  computer  time  demands,  the  Conductive  Media  DC  application  mode  was  not  active  in 
this  model,  i.e.,  only  6  application  modes  were  used.  Adsorption  on  an  initially  clean 
adsorbent  bed  was  simulated.  The  model  was  solved  for  7500  s,  until  complete  saturation  of 
the  adsorbent  bed  was  achieved.  The  results  presented  here  were  obtained  for  inlet 
concentration  0.001  mol/mol  and  inert  gas  flow  rate  0.03  mol/s.  All  other  parameters  used  in 
the  model  are  listed  in  Table  1. 

Velocity  distribution 


Time=7500  Slice:  Velocity  Reid  Arrow:  Velocity  field  Max:  66.442 


The  velocity  during  adsorption  is  practically  not 
changing,  so  in  Figure  10  we  show  the  slice  plot 
corresponding  to  only  one  time,  t=7500  s.  In  Figure  11, 
the  2-D  plots  of  the  velocities,  corresponding  to  planes  A 
and  B  defined  in  Figure  5  are  shown,  and  in  Figure  12,  a 
1-D  plot  showing  the  time  profiles  of  the  velocities  in  5 
points  at  the  intersection  of  planes  A  and  B,  also  defined 
in  Figure  5.  The  dashed  lines  in  Figure  12  correspond  to 
Step  2  of  the  solution  procedure  defined  in  Section  3.3. 
(solution  of  the  model  with  fixed  velocities).  These  lines 
were  shown  in  order  to  check  how  big  the  discrepancies 
between  the  results  obtained  in  Step  2  and  Step  3,  were. 


Figure  10.  Velocity  distribution  at 
the  end  of  adsorption,  in  one  quarter 
of  the  4-cartridge  adsorber 
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t=2000s,  Velocity  wg  [cm/s] 


Max:  47.292 


t=2000s.  Velocity  wg  [cm/s] 


Max:  30.225 


Figure  11.  The  velocity  distributions  during  adsorption  corresponding  to  cross-sections  A  and  B  (Fig.  9)  and 
/=2000  s 


Figure  12.  Time  dependence  of  the  velocity  during  adsorption 
in  5  points  in  the  adsorber  (Fig.  9) 


Gas  concentration  distributions 

Figure  13  shows  the  gas  concentrations  corresponding  to  60,  2000  and  7500  s,  in  3-D  plots. 
Figure  14  shows  the  gas  concentration  distributions  across  the  cross-sections  A  and  B.  These 
results  correspond  to  t=2000  s.  And  finally,  in  Figure  15,  we  show  the  time  profiles  of  the  gas 
phase  concentrations  in  the  five  points  defined  in  Figure  9,  together  with  the  curves  obtained 
with  the  fixed  gas  velocities,  in  Step  2.  Figure  15  shows  that  the  difference  between  the 
results  obtained  in  Step  1  and  Step  2  are  rather  small.  The  blue  lines,  corresponding  to  the 
center  of  the  adsorber  are  not  visible,  as  they  overlap  with  the  violet  lines  corresponding  to 
the  point  x=y=3.29  cm.  This  is  a  result  of  high  diffusion  coefficients  used  in  the  simulations, 
resulting  with  practically  uniform  concentration  in  the  space  around  the  cartridges. 
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Min:  0 


Figure  13.  Concentration  in  the  gas  phase  in  one  quarter  of  the  4-cartridge  adsorber,  during  adsorption 
(from  left  to  right  t=6 0,  2000  and  7500  s) 


Figure  14.  The  gas  concentrations  during  adsorption  corresponding  to  cross-sections  A  and  B  (Fig.  9)  and  £=2000  s 


Figure  15.  Time  dependence  of  the  gas  concentration  during  adsorption 
in  5  points  in  the  adsorber  (Fig.  9) 
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Solid  concentration  distributions 


The  simulated  concentrations  in  the  solid  phase  corresponding  to  60,  2000  and  7500  s,  are 
first  shown  in  3-D  plots,  in  Figure  16.  The  2-D  pictures  of  the  solid  concentrations  at  cross- 
sections  A  and  B,  for  t=2000  s,  are  given  in  Figure  17  and  their  time  profiles  in  Figure  18. 
Naturally,  only  the  results  corresponding  to  the  two  points  inside  the  adsorbent  bed  are  shown 
in  this  figure.  As  previously,  the  results  of  the  approximate  solution  obtained  in  Step  2  are 
also  shown  in  Figure  18,  as  dashed  lines.  The  differences  are  again  rather  small. 


r  me=60  S(oe:  Canae'fc'A.aT,  q  A"ow:  Velocity  f-eld 


T*ne=2000  SCce:  Ganaeifci^at  q  A"ow:  Vetabtf  field 


Tmc=73X  Stae:  CnnaMt.^a-i,  q 


Min:  0 


Figure  16.  Concentrations  in  the  sold  phase  in  one  quarter  of  the  4-cartridge  adsorber,  during  adsorption 
(from  left  to  right  t=6 0,  2000  and  7500  s) 


Figure  17.  The  solid  concentrations  during  adsorption  corresponding  to  cross-sections  A  and  B  (Fig.  9) 
and  t= 2000  s 
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Figure  18.  Time  dependence  of  solid  concentration  during  adsorption 
in  2  points  in  the  adsorbent  bed  (Fig.  9) 


Gas  temperature  distribution 

The  3-D  plots  of  the  temperatures  of  the  gas  phase  corresponding  to  60,  2000  and  7500  s  are 
given  in  Figure  19.  Figure  20  shows  the  gas  temperatures  across  cross-sections  A  and  B,  for 
t= 2000  s,  and  Figure  21  the  time  profiles  of  the  gas  temperatures  in  the  five  points  defined  in 
Figure  5.  The  dashed  lines  again  correspond  to  the  solution  obtained  in  Step  2,  with  fixed 
velocities.  Although  the  differences  between  the  solid  and  the  dashed  lines  seem  remarkable, 
actually  they  are  less  than  2  K.  The  lines  corresponding  to  the  center  of  the  adsorber  and  to 
the  point  closest  to  the  adsorber  wall,  again  overlap  (the  blue  lines  are  not  visible). 


Tnie=®)  Sloe:  Ten  pe>4tuie  Velocity  field 


Min:  290 


Figure  19.  The  temperatures  of  the  gas  phase  in  one  quarter  of  the  4-cartridge  adsorber,  during  adsorption 
(from  left  to  right  t=6 0,  2000  and  7500  s) 
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in  5  points  in  the  adsorber  (Fig.  9) 


Solid  temperature  distribution 

The  simulated  temperatures  of  the  solid  phase  are  shown  in  Figure  22  (the  3-D  plots),  Figure 
23  (the  2-D  plots)  and  Figure  24  (the  time  profiles  in  the  1-D  plot).  The  presentation  is 
completely  analogous  to  the  one  used  for  the  solid  concentrations,  in  Figures  16,  17  and  18. 
The  discrepancies  between  the  final  solution  and  the  solution  obtained  in  Step  2  are  of  the 
same  order  of  magnitude  as  for  the  gas  temperature. 
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Min:  290 


Figure  22.  Temperatures  of  the  solid  phase  in  one  quarter  of  the  4-cartridge  adsorber,  during  adsorption 
(from  left  to  right  t=6 0,  2000  and  7500  s) 


l=2000s,  Temperature  Ts  [K]  Max :  296 .233  t=2000s.  Temperature  Ts  [K]  Max :  296 . 104 


Figure  23.  The  solid  temperature  during  adsorption  corresponding  to  cross-sections  A  and  B  (Fig.  9)  and 
1= 2000  s 


Figure  24.  Time  dependence  of  the  solid  temperature  during  ads.  in  2  points  in  the  adsorbent  bed  (Fig.  9) 
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4.2.  Simulation  of  desorption  without  condensation 

The  desorption  process  before  the  start  of  condensation  at  the  adsorber  wall  was  simulated 
using  Model_D3.  All  7  application  modes  defined  in  Section  2.2.  were  active  in  this  model. 
The  simulation  results  presented  in  this  report  correspond  to  desorption  from  a  previously 
saturated  adsorption  bed,  in  equilibrium  with  the  gas  phase  of  concentration  0.001  mol/mol. 
The  inlet  gas  phase  is  pure  inert,  with  flow  rate  0.01  mol/s,  and  the  voltage  is  3  V.  All  other 
simulation  parameters  are  defined  in  Table  1.  For  the  chosen  set  of  parameters,  at  t=6 0  s  the 
gas  concentration  at  the  adsorber  wall  reaches  0.0147  mol/mol,  which  is  the  saturation 
concentration  of  MEK  corresponding  to  the  wall  temperature  of  293.15  K.  Although  the 
simulation  was  performed  for  longer  times,  only  the  results  for  times  from  0  to  60  s 
correspond  to  this  stage  of  the  TSA  cycle. 

The  simulation  results  for  desorption  without  condensation  are  presented  in  an  analogous  way 
as  for  adsorption.  The  3-D  plots  are  shown  for  t=  10,  40  and  60  s  (the  velocity  only  for  60  s), 
and  the  2-D  plots  for  t= 40  s.  The  results  are  given  in  Figures  25  to  39. 


Velocity  distributions 


Time=60  Slices  Velocity  field  Arrow:  Velocity  field 


Max:  17 .5 


16 


14 


12 


Figure  25.  Velocity  distribution  during  desorption  in  one  quarter  of 
the  4-cartridge  adsorber,  for  t=6 0  s 


t=40s.  Velocity  wg  [cm/s] 


MinsO 


t=40s.  Velocity  wg  [cm/s] 


MinsO 


Figure  26.  Velocity  distribution  during  desorption  corresponding  to  cross-sections  A  and  B  (Fig.  9)  and  t=  40  s 
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Figure  27.  Time  dependence  of  the  velocity  during  adsorption  in  5  points  in  the  adsorber  (Fig.  5) 


Gas  concentration  distribution 


Time=10  Slice!  Concentration,  c  Arrow!  Velocity  field 


Time=40  Slice!  Concentration,  c  Arrow:  Velocity  field 


Time=60  Slice:  Concentration,  c  Arrow:  Velocity  field  Max:  0,0147 


Min:  0 


Figure  28.  Concentration  of  the  gas  phase  in  one  quarter  of  the  4-cartridge  adsorber,  during  desorption 
(from  left  to  right  *=10,  40  and  60  s) 


Figure  29.  The  gas  concentration  during  desorption  at  cross-sections  A  and  B  (Fig.  9)  and  for  *=40  s 
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Figure  30.  Time  dependence  of  the  gas  concentration  during  desorption  in  5  points  in  the  adsorber  (Fig.  9) 


Solid  concentration  distribution 


Min:  3.667e-3 


Figure  31.  Concentration  of  the  solid  phase  in  one  quarter  of  the  4-cartridge  adsorber,  during  desorption 
(from  left  to  right  *=10,  40  and  60  s) 


t=40s,  Concentration  q  [mol/g]  Max:  4.865e-3 

xlO'3 


Figure  32.  The  solid  concentration  during  desorption  corresponding  to  cross-sections  A  and  B  (Fig.  9)  and  t=6 0  s 
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Figure  33.  Time  dependence  of  the  solid  concentration  during  desorption  in  2  points  in  the  adsorbtion  bed  (Fig.  9) 


Gas  temperature  distribution 


Time=10  Slice:  Temperature  Arrow:  Velocity  field 


Time=40  Slice:  Temperature  Arrow:  Velocity  field 


Figure  34.  Temperature  of  the  gas  phase  in  one  quarter  of  the  4-cartridge  adsorber,  during  desorption 
(from  left  to  right  2=10,  40  and  60  s) 


t=40s,  Temperature  Tg  [K]  Max:  340.38  t=40s,  Temperature  Tg  [K]  Max:  340.38 


Figure  35.  The  gas  temperature  during  desorption,  corresponding  to  cross-sections  A  and  B  (Fig.  9)  and  *=40  s 
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Figure  36.  Time  dependence  of  the  gas  temperature  during  desorption  in  5  points  in  the  adsorber  (Fig.  9) 


Solid  temperature  distribution 


Figure  37.  Temperature  of  the  solid  phase  in  one  quarter  of  the  4-cartridge  adsorber,  during  desorption 
(from  left  to  right  t=  10,  40  and  60  s) 


t=40s,  Temperature  Ts  [K]  Max :  341 .23 


t=40s.  Temperature  Ts  [K]  Max:  341.23 


Figure  38.  The  solid  temperature  during  desorption,  corresponding  to  cross-sections  A  and  B  (Fig.  9)  and  t= 40  s 
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Figure  39.  Time  dependence  of  the  solid  temperature  during  desorption,  in  2  points  in  the  adsorbent  bed  (Fig.  9) 


4.3.  Simulation  of  desorption  accompanied  with  condensation 

The  desorption  process  accompanied  with  condensation  of  the  desorbed  vapor  at  the  adsorber 
wall  was  simulated  using  Model_DC3.  All  7  application  modes  defined  in  Section  2.2.  were 
again  active  in  this  model.  The  simulation  of  this  stage  was  performed  as  a  continuation  of  the 
previous  one  -  desorption  without  condensation,  i.e.  the  stored  solution  of  Model_D3  for  t=60 
s  was  used  for  the  initial  values  of  the  variables  for  Model_DC3.  Model_D3  and  Model_DC3 
differ  only  in  the  boundary  conditions  for  the  mass  and  heat  transfer  equations  at  the  adsorber 
wall  (equations  (29)  and  (29a)).  Faster  convergence  was  obtained  if  the  change  of  these 
boundary  conditions  was  performed  in  two  steps:  first  only  for  the  mass  transfer  equation  and 
then  for  the  heat  transfer  equation,  as  well.  All  model  parameters,  including  the  flow  rate  and 
voltage  were  the  same  as  in  the  previous  stage. 

As  mentioned  in  Section  3.3.,  the  solution  of  the  complete  Model_DC3  did  not  converge, 
even  with  increased  tolerances.  Because  of  that,  the  simulation  results  presented  here  were 
obtained  with  fixed  velocities  (the  final  values  of  the  velocities  obtained  with  Model_D3  for 
t=6 0  s  were  used).  The  simulation  results  for  adsorption  and  desorption  without  condensation 
show  that  the  solution  obtained  with  fixed  velocities  don’t  differ  very  much  from  the  final 
solution  of  the  integral  models,  so  the  approximate  solution  obtained  for  fixed  velocities  is 
acceptable. 

The  simulation  results  for  desorption  accompanied  with  condensation  are  presented  in  an 
analogous  way  as  for  adsorption  and  desorption  without  condensation.  The  only  difference  is 
in  the  time  profiles,  as  the  solution  of  the  integral  model  with  variable  velocities  is  not 
available.  Also,  as  the  velocities  are  kept  constant,  the  velocity  distribution  is  the  same  as  for 
desorption  without  condensation  at  t=  60  s.  so  it  is  not  shown  here.  The  3-D  plots  are  shown 
for  t= 60,  140  and  330  s,  and  the  2-D  plots  for  t=140  s.  The  results  are  given  in  Figures  40  to 
51. 
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Gas  concentration  distribution 


T*ne=GO  Sloe:  CanceflUdUa-v  c  A"aw:  Vetaot y  field 


Figure  40.  Temperature  of  the  gas  phase  in  one  quarter  of  the  4-cartridge  adsorber,  during  desorption  with  condensation 
(from  left  to  right  t=6 0,  140  and  330  s) 


<=  140s,  Concentration  c  [mol/mol]  Max:  0.0176  (=140s.  Concentration  c  [mol/mol]  Max:  0.0172 


Figure  41.  The  gas  concentration  during  desorption  with  condensation,  corresponding  to  cross-sections  A  and  B 
(Fig.  9)  and  t=  40  s 


Figure  42.  Time  dependence  of  the  gas  concentration  during  desorption  with  condensation,  in  5  points  in  the  adsorber  (Fig.  9) 
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Solid  concentration  distributions 


Figure  43.  Concentration  of  the  solid  phase  in  one  quarter  of  the  4-cartridge  adsorber,  during  desorption  with  condensation 
(from  left  to  right  t=6 0,  140  and  330  s) 


t=140s,  Concentration  q  [mol/g] 


Mins  3,65e-3 


t=140s,  Concentration  q  [mol/g]  Max!  4.568e-3 

xlO'3 


Figure  44.  The  solid  concentration  during  desorption  with  condensation,  corresponding  to  cross-sections  A  and  B 
(Fig.  9)  and  £=140  s 


Figure  45.  Time  dependence  of  the  solid  concentration  during  desorption  with  condensation,  in  2  points  in  the  adsorbtion  bed 
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Gas  temperature  distributions 


Figure  46.  Temperature  of  the  gas  phase  in  one  quarter  of  the  4-cartridge  adsorber,  during  desorption  with  condensation 
(from  left  to  right  t=6 0,  140  and  330  s) 


Figure  48.  Time  dependence  of  the  gas  temperature  during  desorption  with  condensation  in  2  points  in  the  adsorbtion  bed 
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Solid  temperature  distributions 


Tme=flO  Stce:  Tern  peidtuie  Amwr:  Vetaotf  field 


Figure  49.  Temperature  of  the  solid  phase  in  one  quarter  of  the  4-cartridge  adsorber,  during  desorption  with  condensation 
(from  left  to  right  t=6 0,  140  and  330  s) 


Figure  50.  The  solid  temperature  during  desorption  with  condensation,  corresponding  to  cross-sections  A  and  B 
(Fig.  9)  and  *=40  s 


Figure  51.  Time  dependence  of  the  solid  temperature  during  desorption  with  condensation,  in  2  points  in  the  adsorbtion  bed 
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The  condensation  rate  has  also  been  calculated,  using  equations  (30)  and  (31).  Its  time 
dependence  is  shown  in  Figure  52. 


Figure  52.  Condensation  rate  in  the  4-cartridge  adsorber 


Figure  52  shows  that  the  condensation  starts  at  60  s  and  ends  at  330  s.  Integration  of  the 
condensation  rate  over  time  gives  the  total  amount  of  the  condensed  liquid: 

Lcond  =  1-415  xl0“2  mol 


The  used  electric  power  per  cartridge  was  also  calculated,  as  defined  by  equation  (33)  and 
(34).  It  was  calculated  for  both  stages  of  desorption,  before  and  after  the  start  of  condensation. 
It  is  shown  in  Figure  53. 


Figure  53.  The  used  electric  power  during  desorption,  per  cartridge 


Integration  of  the  electric  power  over  time  gives  the  spent  electric  energy  during  desorption. 
For  the  whole  adsorber  with  four  cartridges,  this  energy  is: 

E  =  110.272  kJ 
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5.  Conclusions 


The  objective  of  this  phase  of  the  project:  to  model  the  TSA  system  with  an  adsorber  with 
four  parallel,  cartridge-type  adsorbers,  raised  a  very  difficult  and  demanding  task,  even  when 
using  an  excellent  modeling  tool,  such  as  COMSOL  Multiphysics.  It  was  necessary  to  build 
3-D  models  of  a  very  complex  system,  with  complex  geometry  and  a  large  number  of 
interacting  phenomena  and  processes.  Even  when  the  problem  was  reduced  to  modeling  of 
only  one  quarter  of  the  adsorber  (assuming  that  all  four  cartridges  in  the  adsorber  were 
identical),  seven  interacting  application  modes  needed  to  be  incorporated  into  a  3-D  model 
geometry  consisting  of  3  sub-domains.  This  resulted  with  a  very  large  number  of  degrees  of 
freedom,  i.e.,  with  a  large  number  of  coupled  nonlinear  differential  equations  that  needed  to 
be  solved  numerically. 

The  first  problem  associated  with  that  was  a  very  high  demand  regarding  computer  memory. 
A  rather  coarse  mesh  had  to  be  used  in  order  to  enable  solving  of  the  models  on  a  personal 
computer  with  4  GB  RAM. 

The  second  problem  was  diverging  of  the  solution.  The  complete  models  did  not  converge. 
For  that  reason,  a  step-by-step  solution  procedure  had  to  be  developed.  For  the  model  of 
desorption  with  condensation,  even  this  procedure  failed,  and  only  an  approximate  solution  of 
the  model  with  fixed  velocities  could  be  obtained.  We  believe  that  this  procedure  could  help 
other  researchers  in  solving  their  difficult  3-D  problems. 

The  third  problem  was  in  relatively  long  computing  times.  For  example,  when  solving  the 
model  for  adsorption,  we  needed  -200  s  for  Step  1,  -1 150  s  for  Step  2  and  -4300  s  for  Step  3. 

Nevertheless,  in  spite  of  all  these  problems,  we  succeeded  in  postulating  and  solving  the 
models  describing  all  three  stages  of  the  TSA  process  in  the  4-cartridge  adsorber. 

We  believe  that  the  convergence  problems  are  partly  associated  with  the  coarse  mesh  used, 
and  that  using  a  finer  mesh  would  improve  the  convergence,  on  account  of  a  larger  computer 
memory  and  longer  computer  times.  Of  course,  this  would  demand  using  of  a  much  faster 
computer  with  more  memory. 

With  this  phase  of  the  project  we  completed  mathematical  modeling  of  the  adsorber 
configurations  developed  in  the  group  of  Prof.  Rood.  These  models  can  be  used  for 
simulation,  optimization  of  the  process  parameters  and  adsorber  geometry,  but  also  for 
comparison  between  different  configurations.  Nevertheless,  for  practical  applications  of  these 
models,  it  would  be  necessary  to  have  accurate  physical  and  transport  parameters  for  the 
particular  system.  Usually,  such  parameters  have  to  be  estimated  from  experimental  data. 
Also,  more  accurate  solution  of  the  mathematical  models,  obtained  with  finer  mesh  and  lower 
tolerances,  would  be  needed,  which  is  directly  related  to  the  computer  hardware.  This  is 
especially  important  for  the  3-D  models  developed  in  this  phase  of  the  project. 

The  set  of  mathematical  models  developed  in  the  three  phases  of  the  project  offer  a  good  basis 
for  optimization  of  the  TSA  processes  with  electrothermal  desorption. 
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Nomenclature 


a  (m2/m3) 
b  (K-1) 

C  (mol/mol) 
C*(mol/mol) 

Csal  (mol/mol) 
cpg  (J/mol/K) 
cpi  (J/mol/K) 
cps  (J/g/K) 
cpv  (J/mol/K) 

Dm  (mol/cm/s) 
Ahg  (W/K/cm) 
Ahs  (W/K/cm) 

E  (J/mol) 
g  (cm/s2) 

G  (mol/s) 

A  (cm) 
hb  (J/cm2/K) 
hs\  (J/cm2/K) 
hs2  (J/cm2/K) 
hws  (J/cm2/K) 

J  (A/cm2) 

Jcond  (mol/cm2/s) 
k  (cm2) 

km  (mol/cm2/s) 
km  i  (mol/cm2/s) 
km 2  (mol/cm2/s) 

4 and  (m°1/s) 

Econd  (mol) 

P  (Pa) 

Pa  (Pa) 

Pc  (Pa) 

4  (Pa) 

QA  W) 

QelO) 
q  (mol/g) 

Rg  (J/mol/K) 


-  Specific  surface  area 

-  Temperature  coefficient  of  the  bed  electrical  resistivity 

-  Adsorbate  concentration  in  the  gas  phase 

-  Adsorbate  concentration  in  the  gas  phase  in  equilibrium  with  the  solid  phase 

-  Saturation  concentration 

-  Specific  heat  capacity  of  the  inert  gas 

-  Specific  heat  capacity  of  liquid  adsorbate 

-  Heat  capacity  of  the  solid  phase 

-  Specific  heat  capacity  of  gaseous  adsorbate 

-  Mass  transfer  dispersion  coefficient 

-  Heat  diffusivity  of  the  gas  phase 

-  Heat  diffusivity  of  the  solid  phase 

-  Adsorption  energy  of  the  adsorbate  (D-R  eq.) 

-  Gravitation  constant 

-  Flow  rate  of  the  inert  gas 

-  Bed  axial  dimension  (Fig.  7) 

-  Solid  to  gas  heat  transfer  coefficient 

-  Heat  transfer  coefficient  from  the  solid  phase  to  the 

-  Heat  transfer  coefficient  from  the  solid  phase  to  the 

-  Gas  to  ambient  heat  transfer  coefficient  (heat  losses) 

-  Current  density 

-  Condensation  flux 

-  Asorbent  bed  permeability 

-  Mass  transfer  coefficient  in  the  adsorbent  bed 

-  Mass  transfer  coefficient  between  the  adsorbent  bed  and  the  gas  in  the  central  tube 

-  Mass  transfer  coefficient  between  the  adsorbent  bed  and  the  gas  in  the  annular  tube 

-  Flow-rate  of  the  condensed  liquid 

-  Total  amount  of  the  condensed  liquid 

-  Gas  pressure 

-  Adsorbate  partial  pressure 

-  Critical  pressure 

-  Adsorbate  saturation  pressure 

-  Rate  of  heat  generation  (equal  to  electric  power) 

-  Heat  generation  (equal  to  electric  energy) 

-  Adsorbate  concentration  in  the  solid  phase 

-  Gas  constant 


gas  phase  in  the  central  tube(s) 
gas  phase  in  the  annular  tube 
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r\  (cm)  -  radius  of  the  central  tube  (Fig.  7) 

r2  (cm)  -  outer  radius  of  the  adsorbent  bed  (Fig.  7) 

rx  (cm)  -  radius  of  the  adsorber  vessel  (Fig.  7) 

Fa  (K)  -  Ambient  temperature 

Fg  (K)  -  Gas  phase  temperature 

Tc  (K)  -  Critical  temperature 

7r  (K)  -  Referent  temperature 

Fs  (K)  -  Solid  phase  temperature 

Tw  (K)  -  Wall  temperature 

t  (s)  -  Time 

U  (V)  -  Electric  potential 

U0  (V)  -  Supply  voltage 

u  (cm/s)  -  x-component  of  the  gas  velocity 

VP  a,  VPB,  VPC ,  VPD  -  Wagner  constants 

v  (cm/s)  -  y-component  of  the  gas  velocity 

W  (cm3/g)  -  Adsorbate  volume  per  1  g  of  the  adsorbent  mass 

Wo  (cm3/g)  -  Total  volume  of  micropores  (D-R  equation) 

w  (cm/s)  -  z-component  of  the  gas  velocity 

x  (cm)  -  x  -  coordinate 

x0  (cm)  -  position  of  the  cartridge  center  in  the  x-direction 

y  (cm)  -y  -  coordinate 

yo  (cm)  -  position  of  the  cartridge  center  in  the  y-direction 

z  (cm)  -  z  -  coordinate 


Greek  letters: 

A//ads  (J/mol) 

-  Heat  of  adsorption 

A//COnd  (J/mol) 

-  Heat  of  condensation 

-  Bed  porosity 

jli  (Pa/s) 

-  Dynamic  viscosity 

p  (Qcm) 

-  Electric  resistivity 

p0  (Qcm) 

-  Electric  resistivity  at  referent  temperature  TR 

pg  (mol/cm3) 

-  Gas  phase  density 

pb(g/cm3) 

-  Adsorbent  bed  density 

pA(g/cm3) 

-  Adsorbate  density 
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Subscripts: 


b 

-  bed 

g 

-  gas 

in 

-  inlet 

i,o 

-  inner  and  outer  tube 

out 

-  outlet 

P 

-  previous  (initial) 

s 

-  solid  phase 

it 

-  inner  (central)  tube 

ot 

-  outer  (annular)  tube 

X 

-  x-direction 

y 

-  ^-direction 

z 

-  z-direction 
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